Different scenarios of topological phase transitions in homogeneous neutron matter 
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We study different scenarios of topological phase transitions in the vicinity of it condensation 
point in neutron matter. The transitions occur between the Fermi liquid state and a topologically 
different one with two sheets of the Fermi surface. Two possibilities of a rearrangement of quasi- 
particle degrees of freedom are shown: the first order topological phase transition and the second 
order one. The order of the phase transition is found to be strongly dependent on the value of the 
critical wave vector of the soft 7r° mode. The thermodynamics of the system is also studied. It is 
shown that the topology of the quasiparticle momentum distribution is mainly determined by the 
neutron matter density, while the temperature T is essential in a narrow density region. A simple 
explanation of the first order topological phase transition at T = is given. 

PACS numbers: 21.65.-f, 26.60.-c, 71. 10. Ay 



I. INTRODUCTION 

A dense neutron matter is an example of the system 
which can exhibit strong correlations of quasiparticles 
(QP) degrees of freedom. The enhancement of QP inter- 
action is revealed 0, in the vicinity of it condensation 
point (PCP) 3 where the system acquires a condensate 
of spin-isospin excitations with quantum numbers of the 
7T° meson. The microscopic calculations [J, |5j predict the 
critical density value p c ~ 0.2 fm -3 for the PCP in neu- 
tron matter. This means that in a typical neutron star 
with the center density of 0.5 — 1.0 fm~ 3 a significant part 
of stellar bulk should exist in a pion condensate phase. 

The non-Fermi liquid topology of the ground state QP 
momentum distribution n(p) in the vicinity of PCP was 
first uncovered in 0. With the increase of the den- 
sity p towards the PCP value p c the QP spectrum e(p) 
(measured from the chemical potential p) becomes a non- 
monotonic function and at certain density pi, manages to 
touch the momentum axis at some point pb (see Fig.Q]left 
panel). This situation is associated with a quantum crit- 
ical point (QCP) 1 at which the single particle density 
of states diverges [7j. Beyond the QCP, the Landau state 
with the Fermi step QP distribution tifl(p) = 0(j>f ~ p) 
becomes unstable, violating the necessary stability con- 
dition 




FIG. 1. Left panel: QCP at the pb density value. Right 
panel: Beyond the QCP. 



A new state appears with a "bubble" (unoccupied re- 
gion) in the momentum distribution, which has several 
sheets of the Fermi surface (see Fig. [T] right panel). To 
our knowledge, such state was first discussed in J8| and 
the one bubble existence problem was studied in |9T[ll|. 
Further developments concerning multi-bubble states can 
be found in [l2|, [H| ■ 

All reconstructions of the QP momentum distribu- 
tion in strongly correlated Fermi systems that cause a 
change of the Fermi surface topology are generally re- 



ferred [14] as topological phase transitions (TPT). It 
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n(p)}) Sn(p) dv > 0. 



(1) 



Here dv = d 3 p/(27r) 3 is the volume element in the mo- 
mentum space and the factor two stands for summation 
over spin projections. The constraint ([T|) requires a posi- 
tive change of the total system energy E[n] for any admis- 
sible variation 8n(p) of the QP momentum distribution 
that conserves the density 



2 J 5n(p) dv = 0. 



(2) 



should be noted that besides the bubble scenario, there is 
another type of TPT which is called the fermion conden- 
sation [15l4l7j . In such TPT the QP spectrum acquires 
a flat band s(p) = 0, p S [Pi,p/], and the Fermi surface 
changes its dimension. The relation between the bubble 
and the fermion condensation scenarios is discussed in 
18]. 

Going back to the neutron matter problem, we stress 
that a bubble formation beyond the QCP is quite im- 
portant for the physics of neutron stars cooling. In- 
deed, a new sheet of the neutron Fermi surface, with a 
low value of the corresponding Fermi momentum, plays 
an important role for an operation of the direct Urea 
processes: n — > p + e~ + i7 e , p + e~ n + v e - In 
cores of neutron stars this processes are generally con- 
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sidered to be forbidden due to the kinematic restric- 
tion for the Fermi momenta of the involving particles: 
pp n < Pf p + Pf c ■ In the typical density range of the core 
~ 1 — 2 po (po — 0.16 fm~ 3 is the normal nuclear density) 
the proton fraction not exceeds 6 — 8%, and the r.h.s. of 
the kinematic equality is estimated by ~ 0.8 pF n - The 
appearance of the new sheet of the Fermi surface at the 
point py < 0.8pp n allows to agree the kinematic re- 
striction. This mechanism was considered in @, H(| as a 
possibility of the enhanced cooling of some neutron stars 
(e.g. Vela, Geminga, 3C58). 



II. QUASIPARTICLE APPROACH NEAR PION 
CONDENSATION POINT 

A. General relations 

We begin with a concise discussion of a method for 
describing neutron matter near the PCP. It is based on 
an implementation of the Landau-Migdal QP approach 
to strongly correlated Fermi systems that is reviewed in 
details 21]. The QP spectrum and the momentum dis- 
tribution can be found from the set of equations: 



m 



0s(p) 
dp 

n(p,T) = [l 
u(p) dv = p . 



, dn(p') , 
f(p,p) „ . dv . 



dp' 



(3) 
(4) 
(5) 



The first equation is the Landau relation for the QP spec- 
trum and the interaction amplitude / in homogeneous 
systems. It is based on the translational invariance [22| . 
The quantity m is the free neutron mass. The second 
equation is the Fermi-Dirac statistic formula, where T 
is the temperature. The last one is the normalization 
condition including the neutron matter density p. 

The QP interaction amplitude is identified [22J with 
the w-limit of the vertex function T 

/(Pi,o-i;p2,<7 2 ) = Z 2 T^ /3rtS (p 1 ,p 2 ) 

= Z 2 ■ lim r Q( g ;75 (pi,p 2 ,k,w) , (6) 

where Z is the residue of the single particle Green func- 
tion and <Tq 7 stands for Pauli spin matrices. Accord- 
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FIG. 2. The main contribution to the QP interaction 
amplitude. 



ingly to [l|, the most singular contribution to the ver- 
tex function near the PCP comes from an exchange by 
a soft spin-isospin collective mode. The direct and the 
exchange graphs of the interaction are shown in Fig. [5] 
After taking the k — > limit, the contribution of the 
exchange graph (the second one) still shows a strong de- 
pendence on the relative momentum pi — P2- Thus, the 
QP interaction amplitude reads [2j: 

/(pi,eri;p2,<x 2 ) ~ (JasDJp^fa -p 2 ,w = 0-p) , (7) 

where J a s is the interaction vertex of nucleons and pions 
in neutron media, D is the ir° propagator, and the general 
argument of operators is out of brackets. 



B. Topological phase transitions 

As it is discussed in ||, the spectrum u)(q) of 7r°-like 
QP excitations in neutron matter is given by a particular 
poles branch of the ir° propagator. The behavior of the 
branch depends on the density p as a parameter. Pion 
condensation occurs at the p c value when the excitation 
energy vanishes, uj(q c ) — 0, together with its derivative, 
du)(q c )/dq = 0, at a certain momentum q c . As a conse- 
quence, at the PCP the following conditions are valid: 

(g c ,0;p c )=0, — 



D~ 



0. 



(8) 



Using the Taylor expansion of D^ 1 in the vicinity of 
PCP, the interaction amplitude /, entering in Eq. (|3|), 
can be written [H, [2(| in the form 



/(Pi,P2) = -Sp . 1 Sp . 3 /(pi. 



0-i;p2,o- 2 J 
9 



2 (P) + ((P1-P2) 2 /9 C 2 -1) Z 



(9) 



Here g is an effective coupling constant and n 2 (p) c>c 
(p c — p) is an effective radius in momentum space. Previ- 
ous investigations [23( were focused on the ground state 
topology, while g,q c ,n quantities were regarded phe- 
nomenologically as free parameters. For convenience of 
readers, we present in Fig. [3] a topological phase diagram 
in q C: k variables obtained in that work. The label FL 
corresponds to the Fermi liquid state and LB^, to states 
with i sheets of the Fermi surface. A change of the den- 
sity p leads to a change of the system position (q c , k) on 
the diagram. Transitions between different regions of the 
phase diagram represent TPT that can happen in neu- 
tron matter. Until recently, all such TPT were consid- 
ered [la . |21| to occur continuously according to violation 
of Eq. (fT]). However, it's not the only possibility and an 
attentive investigation [iij] within the model ([9]) showed 
the existence of first order TPT. Such scenario of the 
Fermi surface reconstruction in a homogeneous isotropic 
Fermi system was first revealed (25[ in a model of strongly 
correlated 2D electron gas beyond the QCP. Below, we 
present a detailed analysis of possible TPT for the neu- 
tron matter problem. 
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FIG. 3. The ground state neutron matter phase dia- 
gram near the PCP 0. 
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where = (Sn pj is the neutron Fermi momentum. 
The pion propagator is connected with the polarization 
operator: D' 1 {q,uj) = uj 2 - q 2 - to 2 - U(q,u). The 
microscopic description of the pion polarization operator 
is a quite subtle matter 26j, and we use here the semi- 
microscopical representation [28j 



Tl(q,0) = -q 2 



Xnn(q) 



p 



The first term describes processes of particle-hole excita- 
tions where nucleon-nucleon correlations are taken into 
account by the denominator. The second term is a phe- 
nomenological one and describes A-resonance-hole exci- 
tations. The ^-scattering processes in neutron matter 
are neglected. Finally, one arrives at the formula 



III. SEMI-MICROSCOPIC QP INTERACTION 
AMPLITUDE 



The above discussion of TPT relied on the phenomeno- 
logical description ((9|) of the QP interaction amplitude. 
It is possible to convey neutron matter physics near the 
PCP in a more direct way. A semi-microscopical expres- 
sion for the QP amplitude can be derived by use of the 
microscopic formula (|7|). 

The bare 7r°-nucleon interaction vertex is given (26| by 



7° 

TY t 



if 



NN 



(o-q)r 3 



(10) 



where / ~ 1 is the dimensionless 7r°-nucleon coupling 
constant, m n is the pion mass and T3 is the diagonal 
isospin matrix. The vertex renormalization in neutron 
matter is due to nucleon-nucleon correlations that can 
be described by means of the Landau-Migdal amplitude 
271 



m po 



7T- 



F=F + F't 1 T2 + (G + G't 1 T 2 )<T1<T2. (11) 



1 /3 

Here po = (l.57r 2 /?o) is the Fermi momentum in equi- 
librium nuclear matter and m* is the effective nucleon 
mass. Accordingly to Q, the renormalized static vertex 
reads 



if (o"a/3q)T3 

1 + 9cXnn(i) 



m*p f 



.</. = - r —g(G + G'). 

(12) 

The function Xnn^) * s * ne static susceptibility of free 
QP proportional to the Linhard function: 



X s nn(9) 



f 2 m*p F 



pf / q 



2 2<z \4p 2 F 



1 - q/2p F 



1 + q/2p F 



(13) 



f(q) = 2 Sp tri Sp . 2 /(pi,o-i;p 2 ,<T 2 ) 
/ q 



tt 1 + 9cX S nn (<?) / ml + q 2 + U(q, 0) 



(15) 



where the polarization operator is given by (I14[) . It is 
worth to remark that the momentum dependence is en- 
tirely given through the relative momentum value q — 
|Pi - P2|. 
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FIG. 4. The interaction amplitude (|15[l multiplied by 
JVo = rnpF /f 2 as a function of q/pF- The parametriza- 
tion (|16|l corresponds to curves peaked at q c ~ 0.74jjf, 
while (TiT)! , at q c ~ pf ■ 



We now discuss the parameters entering the amplitude 
f(q). The constants: / = 1, m* = m, G + G' = 1 are 
fixed all along the further discussion, while we use two 
sets of parametrizations for the phenomenological part of 
the polarization operator. The first one corresponds to 
the choice in (28[, where 



^= 0.59 77-4, qA = 2.08m„ 



(16) 



The second parametrization is suggested in this work: 

Pa = 0.97 m\ , q A =i.lm w . (17) 
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Each parameterization reproduces the value p c ~ 0.2 fm -3 
that agrees with the critical density of ir° condensation 
obtained in [3, [f| . The difference is in the corresponding 
value of the critical wave vector q Cl which is not known 
well from microscopic calculations. We also note that p c 
and q c are quite sensitive to the tuning of pa and qA 
parameters. The behavior of the QP interaction ampli- 
tude is shown in Fig. |U The parametrization (fTB| leads 
to q c ~ 0.74pi? (that is less then pp), while (fT7|) . to 
q c ~ pp. Fig. 0] also demonstrates the amplification of 
the amplitude with an increase of the pp value (or the 
density p — p F /3n 2 ). The real divergence, that is ex- 
actly at the PCP, is reached at the point pp ~ 2.602 m T 
and pf ~ 2.598 m, correspondingly for the first and the 
second parametrizations. 



IV. DIFFERENT SCENARIOS OF 
TOPOLOGICAL PHASE TRANSITIONS 

Here we present the results for topological rearrange- 
ments of QP degrees of freedom based on the semi- 
microscopical QP interaction amplitude. The associated 
TPT are considered with respect to the density and the 
temperature changes and turns out to be of first or second 
order depending on 7r° polarization operator parameters 
(finally on the q c value). 

The QP spectrum and the momentum distribution are 
determined by the set of equations (J3]) , d4j) , d5j . The first 
equation permits an integration by parts, due to the as- 
sumed dependence (|15|) on the relative momentum q of 
the interaction amplitude 



:(p) = f^ + J.f(p-p')n(p')dv'. 



(18) 



Then, it is convenient to apply a contracting iterations 
method for the solution of this set of equations. The 
case of T — (when it is addressed) was modeled by 
T = 10- 5 e° F , where e% = p%/2m. 

First, we consider the results obtained within the 
model (IT5|) where q c = 0.74p^ . The corresponding QP 
degrees of freedom evolution, driven by an increase of 
Pp, is displayed in Fig. [5] It has a continuous behav- 
ior according to the QCP scenario of the Fermi surface 
reconstruction. At the critical value p b F — 2.562 m w a 
bifurcation occurs, and a new root of the QP spectrum 
appears at the point pb = 0. Beyond the QCP, the QP 
momentum distribution possesses two sheets of the Fermi 
surface with coordinates p\ and p2 — Pf- The bubble re- 
gion size (which is equal to p\) increases continuously 
from the zero value with the further increase of pf ■ 

Regarding Eq. (ITSl , it is easy to write down the corre- 
sponding energy functional 



E[n}=2 J L- n (p)dv + J f(p-p')n(p)n(jpT) 



dv dv' . 
(19) 
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FIG. 5. The QP spectrum (in units e° F = p%/2m) and 
the momentum distribution evolution with the increase 
of pf within the model (|16[) . The temperature T — 0. 
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FIG. 6. Energy per particle at T = as a function of 
p F . Left panel corresponds to the model flB) , while 
right one, to the model (|17[) . 



According to the continuous picture of the ground state 
evolution (Fig. [5]) , the energy shows a monotonic behav- 
ior as it is plotted on the left panel of Fig. [5] The topo- 
logical transition between the FL and the LB2 states is 
of the second order. 

Now, we turn to the model (fl~6|) where q c ~ pf- 
Let us first consider the right panel of Fig. [SJ It is 
seen that there are two different states in the interval 
2.57 m-K < pf < 2.58 m w . We note that in order to estab- 
lish the coexistence of several solutions of the equations 
(|4"]).([5"],.([T8]) it is steadily necessary [H, [H[ to consider 
different initial conditions for the iteration procedure. 
Fig. [7] represents QP spectra and momentum distribu- 
tions of coexisting states. At the value pp = 2.572 
the ground state is the FL one (panel (a)), while there 
is also a LP>2 state (panel (b)) with a slightly higher en- 
ergy value (see the right panel of Fig. |6|). It is worth to 
emphasize that the bubble region has finite, not negligi- 
ble size p\ ~ 0.6pF- A first order phase transition takes 
place at the pf — 2.573 value. Beyond the transition 
point, the LB 2 state (panel (d)) becomes energetically 
favored over the FL state (panel (c)). Finally, the local 
energy minimum of the functional (119[) , corresponding to 
the FL state, becomes unstable |29( (for pp > 2.581 m„) 
and only the LB2 ground state remains (the right panel 

of Fig. ED . 
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FIG. 7. Panels (a), (b): QP spectra and momentum 
distributions of coexisting states in advance of the first 
order TPT. Panels (c), (d): the same, beyond the TPT 
point. The temperature T = 0. 



Thermodynamic functions of neutron matter with the 
Fermi momentum pp = 2.572 are given by Fig. [5] 
The panel (a) demonstrates that the FL state is lower 
in the energy then the LB 2 one up to the temperature 
T ~ 0.1 ~ 5MeV. At the same time, the panel (b) 
shows that the entropy value 

S[n] = -2 J n(p) logn(p) + (1 - n(p)) log(l - n(p)) , 

(20) 
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FIG. 8. The energy (panel (a)), the entropy (panel (b)), 
the free energy (panel (c)) and the density of states 
(panel (d)) as functions of the temperature for the FL 
and the LB2 states. 
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FIG. 9. QP spectra and momentum distributions of the 
FL and the LB2 states at the point (the black dot A in 
Fig. II 0I below) of the temperature driven TPT. 



corresponding to the LB2 state, grows more rapidly as 
the temperature increases, than the one of the FL state. 
An interplay of the energy and the entropy balance leads 
to a first order TPT that is driven by the temperature. 
The free energy F[n] — E[n] — TS[n] behavior is given 
by the panel (c) and demonstrates the phase transition 



at the temperature T\ ~ 6.2 ■ 10" 
(d) shows the density of states 



-2^0 



e F . Finally, the panel 



Af[r 



dn(p) 

de 



dv 



1 

T 



n(p)(l-n(p))dv. (21) 



To elucidate the temperature dependence, we present 
in Fig.[9]QP spectra and momentum distributions of the 
FL state and the LB2 one at the point T = 6.2 • 10 _2 e^ 
near the phase transition. Thus, the higher values of the 
entropy and the density of states, corresponding to the 
LB2 state, are caused by "melting" of two sheets of the 
Fermi surface, while the FL state has only one sheet. 

The temperature investigations of the QP states be- 
havior, in the region near the point pp = 2.572 (see 
the right panel of Fig. [6]), are gathered by a (pp, T)-phasc 




2,573 



FIG. 10. Phase diagram of neutron matter near the 
line Ti(p F ) of the fist order TPT. The first label cor- 
responds to the thermodynamically favored phase, the 
second one, to a metastable phase. The black dot A 
refers to the QP states shown in Fig. [9] 
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diagram of neutron matter shown in Fig. 1101 The dia- 
gram consists of two regions: the first one corresponds to 
the FL phase of the system, while the other one, to the 
LB2 phase. The regions are separated by the line T\(jpp) 
of the first order TPT. The diagram shows that neutron 
matter state topology is determined by the temperature 
in a quite narrow density interval. 

V. ENERGETICS OF LB 2 STATES 

In this section, we elucidate why the system appears 
to be in the LB 2 phase. The energy functional (fT9|) can 
be rewritten in the form 

E[n] = 2 J |-n(p) dv+^f /(q)S(q; [n]) dv . (22) 

Here the interaction energy part is given by use of the 
structure function 

S(q;[n]) = - / n(p + q)n(p)du. (23) 
P J 

Then, we address below the case of T = 0. The set of 
all possible QP momentum distributions with one or two 
sheets of the Fermi surface (that is the FL or LB2 states) 
is specified by 

MP) = *(P2 - p) - 9(pi - p) , (24) 

pl-p\=P%. (25) 

The last equation follows from the normalization con- 
dition ([5]). This class is referred below as n 2 type QP 
distributions. In this case, the structure function is eval- 
uated explicitly, the corresponding formulas are given in 
the Appendix. It is sufficient, due to (|2"5j) . to deal with 
one parameter. The convenient choice is rj = p 2 — pi 
(the width of the occupied region in the QP momentum 
distribution) that defines distribution parameters 



Pi, 2 = 




The value rj = pp corresponds to pi = 0, p 2 — pf (the 
Fermi step). The decrease of rj leads to a monotonic 
increase of pi,2- 

The behavior of the structure function S(q; [712]) = 
S(q; rj) is shown by Fig. [Tl] (the left panel). 

Its explicit form corresponding to the Landau state 
(rj = pp), is well known 

S(q;PF) = S FL {q) 

=\{ i -£)\ 2+ i-p) e ^-^- w 

It is remarkable that in the other important case of thin 
(■q <C Pf) n-i type QP distributions that can appear by a 
first order TPT, one obtains 

S{q\r,)\ 7l ^ = Tn- (28) 




FIG. 11. Left panel: The structure function of the ri2 
type QP momentum distributions. Several values of 
the wave vector q axe regarded. Right panel: The to- 
tal energy per particle for the thin n 2 type momentum 
distributions. 



This result is demonstrated well by Fig. [TTJ and reflects 
the natural decrease of the momentum distributions over- 
lap in the integral of Eq. (|2"3")l . As a consequence, this 
leads to a reduction of the interaction energy value 

The kinetic energy, on the contrary, shows enhancement 
that is confirmed by the explicit result 

3 

E kin [n 2 ] = 3—5- {p\ — p\) e%p 
^Pf 

The total energy value of 712 type distributions in the 
limit rj <C Pf is given by 

E[n 2 }c(^- + ^js%p, u = U/e FP . (31) 

This function is plotted on Fig. [TT] (the right panel) and 
reveals an evident minimum, that is 

E LB2 = J^ue%p (32) 

at the width value r] c = pp v/2/ (3u) . Finally, we note 
that the -Elb 2 energy value appears to be lower the Fermi 
liquid one 

E FL = g4p +\J f(q)S F L(q) dv , (33) 

when the QP interaction is sufficiently strong. Indeed, 
if one characterizes the QP interaction amplitude by an 
effective coupling constant f(q) oc g, equations (|32p and 
([33| show Blb 2 \fd an d -Efl c>c g behavior at the large 
g limit. 

Thus, the appearance of LB 2 states is explained by the 
interplay between the kinetic and the interaction energy 
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contributions. The existence of this energy minimum, 
generally, does not dependent on whether the Landau 
state is present or not. We remark that this simple ex- 
planation is limited by the 712 set of QP momentum distri- 
butions. A generalization to a more complete class n Q >2 
could in principle reveal an instability of a LB2 state with 
respect to some energetically favored one 24|. The solid 
proof of the existence of the LB2 ground state comes from 
the direct solution of the equations (|4jl , ([5]) , (TT8t . 



VI. CONCLUSION 
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VIII. APPENDIX 



In this article we have considered two scenarios of topo- 
logical phase transitions in homogeneous neutron mat- 
ter. The transitions occur between the Fermi liquid state 
and the other one with two sheets of the Fermi surface. 
The investigation was performed with the use of a semi- 
microscopical expression for the quasiparticle interaction 
amplitude in the vicinity of ir° condensation point. The 
order of the phase transition is confirmed to be strongly 
dependent on the critical wave vector q c in agreement 
with more simplified formulations (24|. The first possi- 
bility for a rearrangement of the quasiparticle degrees of 
freedom is the second order topological phase transition. 
It occurs when q c < pp^u-d it corresponds to a quantum 
critical point scenario @, 0, EH of the Fermi surface re- 
construction. The second possibility studied in this work 
is the first order topological phase transition. This case 
occurs when q c > pf and it is connected with a sudden 
change in the quasiparticle momentum distribution and 
spectrum. The first order topological phase transition 
can be driven either by the density or the temperature 
changes. Thermodynamic functions and the phase dia- 
gram are calculated. It is shown that the influence of the 
temperature on the Fermi surface topology is essential in 
a quite narrow density region. A simple explanation of 
the origin of the first order topological phase transition 
at T = is given. 



The structure function within the 112 set of QP mo- 
mentum distributions reads 

S(q; M) = - / n(p + q)n(p) dv = S n + S22 - 2S 12 . 
P J 

(34) 

Here (4tt/3)p f Sij is a volume of two spheres intersection 
with radii p min = pi and p max = pj, p min < p max , while 
the distance between the centers is equal the q vector 
length. The value of the volume is specified by 

Pf Sij{q;pi,Pj) = {C(q;Pi,Pj) + C{q;Pj,Pi)) 

x 9(p t + Pj - q)9{q+Pi-pj) +PiO{p 3 -q-pi). (35) 

The first term corresponds to an intersection case, while 
the second one, to a complete enclosure of the small 
sphere into the bigger one. The function £ has the form 



h \ 1 ( p 2 + i 2 -p 2 \ 
C(q;Pi,Pj) = j I Pi Y q ) 
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